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I \' ABSTRACT 

^ . 

We have discovered two dusty intervening Mg 11 absorption systems at z 1.3 in 

Q , the Sloan Digital Sky Survey (SDSS) database. The overall spectra of both QSOs are 

5h ' red (u-K>4.5 mag) and are well modelled by the composite QSO spectrum reddened 

[/3 I by the extinction curve from the Large Magellanic Cloud(LMC2) Supershell redshifted 

. to the rest-frame of the Mg 11 systems. In particular, we detect clearly the presence 

of the UV extinction bump at Arost ^ 2175 A. Absorption lines of weak transitions 

^ ; like Si 11AI8O8, Cr iiA2056, Cr ii+Zn iiA2062, Mn iiA2594, Ca iiA3934 and Ti iiA1910 

^ I from these systems are detected even in the low signal-to-noise ratio and low resolution 

' SDSS spectra, suggesting high column densities of these species. The depletion pattern 

T— I , inferred from these absorption lines is consistent with that seen in the cold neutral 

G> ■ medium of the LMC. Using the LMC Ay vs. N(R i) relationship we derive A^(H i)'^ 

' 6 X 10^^ cm~^ in both systems. Metallicities are close to solar. Giant Metrewave Radio 

Q>^ . Telescope (GMRT) observations of these two relatively weak radio loud QSOs ~ 

• 50 mJy) resulted in the detection of 21-cm absorption in both cases. We show that the 

00 I spin temperature of the gas is of the order of or smaller than 500 K. These systems 

. provide a unique opportunity to search for molecules and diffuse interstellar bands at 

L • ■ z>l. 
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1 INTRODUCTION 

Studying the physical conditions in the interstellar medium 
(ISM) at high redshift and the processes that maintain these 
conditions is important for our understanding of how galax- 
ies form and evolve. The presence of dust influences the 
physical state of the gas through photo-electric heating, 
UV shielding, and formation of molecules on the surface of 
grains. However, we know very little about dust properties 
in the ISM at high redshifts. Properties of the dust can be 
derived from extinction curves observed in different astro- 
physical objects. Recently, Noll et al. (2007) found evidence 
for the presence of a moderate UV bump in at least part of 
the population of massive galaxies at z > 1. Studies of dust 
extinction in the circum-burst environment of Gamma Ray 
Bursts indicate that in most cases, a LMC-like extinction 
curve is preferred (e.g. Heng et al. 2008). 
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The depletion of Cr with respect to Zn in intervening 
damped Lyman-a systems (DLAs) shows that dust is indeed 
an important component of the high density gas (Pettini, 
Smith & Hunstead 1994) . A correlation is observed in DLAs 
between metallicity and dust-depletion (Ledoux, Petitjean & 
Srianand 2003) which is confirmed by the higher detection 
rate of II2 in DLAs with higher metallicities (Petitjean et 
al. 2006; Noterdaeme et al. 2008). The corresponding gas is 
cold (T~150 K; Srianand et al. 2005) as expected from mul- 
tiphase ISM models in which the gas with high metallicity 
and dust content has lower kinetic temperature than the gas 
with lower metallicity and dust content (Wolfire et al. 1995). 
However, even in the highest metallicity DLAs typical dust 
signatures like high extinction (i.e 0.16< E(B— V)<0.40), 
2175 A absorption bump or the diffuse interstellar bands 
(DIBs) are not seen. 

Signatures of dust are seen from few intermediate and 
low redshift absorption systems. Both DIBs and 2175A ab- 
sorption bump have been detected in the Zabs ~ 0.524 system 
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Table 1. GMRT observation log and results 



J0BS04e.eH-616911.7 



Source name 


Date 


Time 
(hr) 


Peak Flux" 


rms** 


J0850+5159 


2007 Nov 06 


7.3 


64.0 


1.2 


J0852+3435 


2007 Nov 05 


7.5 


51.2 


2.4 




2007 Nov 30 


6.9 


52.2 


1.3 




2008 Mar 08 


6.2 


50.1 


1.3 



" in units of mJY beam ^ 

^ spectral rms in units of mjy beam~^ channel" 



toward AO 0235+164 (Wolfe & Wills 1977; York et al. 2006; 
Junkkarinen et al. 2004; Kulkarni et al. 2007). This system 
also shows strong 21-cm absorption (Wolfe & Wills 1977) 
and has iV(H l) = 5x10^^ cm^^ E(B-V) = 0.23 (Junkkari- 
nen et al. 2004) and Ts = 220±60 K. Wang et al. (2004) 
have reported the detection of the 2175A absorption bump 
in 3 intermediate redshift [z ^ 1.3) Mg ll systems. Recently, 
Ellison et al. (2008) have detected DIBs in the 2abs= 0.1556 
Ca II systems towards SDSS J001342.44-002412.6. Compos- 
ite spectra have been obtained for different samples of ab- 
sorbers from the Sloan Digital Sky Survey (SDSS). They 
show in a statistical way that dust is present in strong Mg ll 
(York et al. 2006a) and Ca ii systems (Wild, Hewett & Pet- 
tini 2006). In the former case the mean extinction curve is 
similar to the SMC curve with a rising ultraviolet extinction 
below 2175 A with E(B-V)<0.08 and with no evidence of 
an UV bump. In the latter case, evidence for the UV bump is 
marginal and LMC extinction curve provides E(B-V)<0.103 
for different sub-samples. 

An efficient way to reveal cold and dusty gas is to search 
for 21-cm absorption as shown by the high detection rate 
of 21-cm absorption towards red QSOs (Carilli et al. 1998; 
Ishwara- Chandra, Dwarakanath & Anantharamaiah 2003; 
Curran et al. 2006). Multiple lines of sight toward lensed 
QSOs passing through regions producing high extinction 
also show 21-cm and molecular absorption lines (see Wiklind 
& Combes 1996). 

We have recently completed a systematic search for 21- 
cm absorption in a complete sample of 38 Mg ll systems, 
drawn from the SDSS DR5, with redshifts in the range 
1.10 < z < 1.45 corresponding to the frequency coverage 
of 610 MHz feed at GMRT (see Gupta et al. 2007 for early 
results). Using an automatic procedure we have identified 5 
systems in this sample showing strong absorption lines from 
Si II, Zn II, Cr ll, Fe ll, Mg ll and Mg I in front of radio- 
loud QSOs with flux density > 50 mjy. Here, we report 
the detection of 21-cm absorption from two of these systems 
at 2abs~1.3 towards red (u-K>4.5 mag) QSOs. The optical 
spectra of these two quasars possess 2175A dust absorption 
features. 



2 OBSERVATIONS AND ANALYSIS: 

We observed J085042. 21-^515911. 7 (called J0850-H5159) and 
J085244.74+343540.5 (called J0852-H3435) with the GMRT 
to search for 21-cm absorption associated with the Mg ll sys- 
tems at 2abs= 1.3265 and ^abs^ 1.3095 respectively. Log of 
our observations is presented in Table [1] We used 1 MHz 




610.2 

J0asg44.7« * 343640.6 




6I4.B Slfi 615.2 616.4 

Heliocentric frequency (MHz) 

Figure 1. GMRT spectra of J085042. 21+515911. 7 (top panel) 
and J085244. 74+343540. 5 (bottom panel). H I 21-cm absorption 
is detected at 2:^133 = 1.3265 and 1.3095 (see Table 2), respectively. 
The solid lines represent the fits to the overall profiles. Individual 
Gaussian components are shown with dotted lines. Residuals, on 
a scale arbitrarily shifted for clarity, are also plotted. 



Table 2. Details of multiple Gaussian fits to 21-cm absorption 
lines. 



QSO 


Zabs 


At;'' 




N{n 1) ^ 


J0850+5159 
J0852+3435 


1.32674 
1.32692 
1.30919 
1.30945 


49+16 
24+4 
23+7 

63+12 


0.117+0.036 
0.347+0.108 
0.078+0.025 
0.079+0.009 


1.11+0.50 
1.61+0.57 
0.35+0.15 
0.96+0.21 


" FWHM in km 


s-l; in 


units of 10 


19 cm-2 K-i 





baseband bandwidth split into 128 frequency channels (cor- 
responding to a velocity resolution of ~3.8 km s~^) cen- 
tered around the redshifted 21-cm frequency. We observed 
standard flux density calibrators 3C 147 and 3C 286 ev- 
ery 2-3 hours to correct for amplitude and bandpass vari- 
ations. Phase cahbrators (J0834+555 for J0850+5159 and 
J0741+312 for J0852+3435) were also observed approxi- 
mately every 40 min. Data were acquired in both the circular 
polarization channels, RR and LL and reduced in the stan- 
dard way using the Astronomical Image Processing System 
(AIPS) as described in Gupta et al. (2006). GMRT spectra of 
both sources show detectable 21-cm absorption (see Fig. [T|. 
The broad absorption produced by the 2abs= 1.3095 system 
towards SDSS J0852+3435 is conflrmed by repeated GMRT 
observations obtained at three different epochs spread over 4 
months. We have used the combined spectrum for the anal- 
ysis presented here. It is usual procedure to decompose the 
absorption profile into multiple Gaussian components. For a 
Gaussian profile the H I column density and the peak optical 
depth (rp) are related by. 
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Table 3. Rest equivalent widths of UV absorption lines. 



Species Rest equivalent widths (A) 

J0850+5159 J0852+3436 Ca ii 
■^abs =1-3265 Zabs =1-3095 composite 



Ca IIA3934 




1.05 


0, 


,49- 


-1, 


,01 


Si IIA1808 


0.96 


1.06 










Mg nA2798 


4.62 


2.89 


2, 


.19- 


-2, 


,47 


Mg IIA2803 


4.15 


2.80 


2, 


.07- 


-2, 


,23 


Mg IA2853 


2.08 


1.11 


0. 


.72- 


-0, 


,93 


Mn nA2576 


0.91 


1.25 


0, 


.20- 


-0, 


,40 


Mn IIA2594 


0.43 


1.30 


0, 


,15- 


-0, 


,23 


Mn IIA2606 


0.50 


0.63 


0, 


,10- 


-0, 


,12 


Fe IIA2249 


0.58 


< 0.56 


0, 


,09- 


-0, 


,10 


Fe IIA2260 


0.55 


< 0.56 


0, 


,08- 


-0, 


,08 


Fe IIA2344 


1.61 


1.46 


1, 


,36- 


-1, 


,39 


Fe IIA2374 


1.72 


1.47 


0, 


,89- 


-1, 


,01 


Fe IIA2382 


2.19 


2.36 


1, 


,55- 


-1, 


,61 


Fe IIA2586 


1.69 


1.60 


1, 


,27- 


-1, 


,35 


Fe IIA2600 


2.27 


2.10 


1, 


,60- 


-1, 


,80 


Zn Il+Mg IA2026 


0.98 


0.85 


0, 


,12- 


-0, 


,27 


Cr IIA2056 


0.67 


< 0.30 


0, 


,09- 


-0, 


,11 


Cr Il+Zn IIA2062 


1.12 


1.20 


0, 


,11- 


-0, 


,20 


Cr IIA2066 


0.50 


0.39 


0, 


,08- 


-0, 


,08 


Ti nA1910 


< 0.10 


< 0.20 










Ti IIA3242 


0.33 


0.60 


0, 


,09- 


-0, 


,12 


Ti nA3385 


0.63 


0.71 




0.09 





Table 4. Results of BED fitting to the SDSS spectrum. 



Object 




Dust 
model 






J0850+5159 


1.89 


MW 


0.83(0.06)1 


2.75 






SMC 


0.59(0.05)1 
0.51(0.04)2 
0.51(0.04)3 


1.42 
1.60 
1.36 






LMC2 


0.83(0.06)1 
0.73(0.05)2 
0.70(0.05)3 


1.16 
1.43 
1.27 


J0852+3435 


1.65 


MW 


1.20(0.08)1 


2.49 






SMC 


0.75(0.06)1 
0.64(0.05)2 
0.68(0.06)3 


2.09 
2.14 
1.99 






LMC2 


1.10(0.08)1 
0.95(0.08)2 
0.97(0.08)3 


1.42 
1.64 
1.53 



1 using the SDSS QSO composite spectrum (Vanden Berk et al. 
2001); 2 using the HST QSO composite spectrum (Zheng et al. 
1997); 3 using LBQS QSO composite spectrum (Francis et al. 
1991). 

Ar(HI)/c = 1.93 X 10^** Tp Ts Au cm"^ (1) 

Here, Ai; and /c are the FWHM (in km s"^) of the fitted 
Gaussian component, and the covering factor of the absorb- 
ing cloud respectively. Results of the Gaussian fits to the 
21-cm absorption lines are given in Table [2] 

Optical spectra of the QSOs were downloaded from the 
SDSS archive. In Table |3] we summarize the rest equivalent 
widths of various metal absorption lines in the two absorp- 
tion systems and compare them with the range of equivalent 



widths observed for the Ca ll systems by Wild, Hewett & 
Pettini (2006). A glance at this table suggests that the two 
systems discussed here have much stronger absorption lines 
than those associated with the Ca ll systems. 

Next we compute the extinction due to the Mg ll system 
by assuming the reddening of the quasar to be a consequence 
of the presence of dust at Zabs- The optical depth at an 
observed wavelength A can be written as, 

r(A) = 0.92ylA, = OmAvAi^a). (2) 

Here, and Av^ are the extinction in magnitude at 

Aa = A/(l -I- «abs) and at the rest V band of the absorber 
respectively. C(-^a) is the relative extinction at Aa to that 
in the rest V band. We consider ^(A) for the SMC, the 
LMC2 Supershell and the Galaxy (Misselt, Clayton & Gor- 
don 1999; Gordon et al. 2003) in order of increasing UV 
bump strength. 

We have developed a robust minimization code using 
the IDL routine MPFIlQ that uses Levenberg-Marquardt 
technique to get best fit values of Av.^ and normalization 
of the flux scale compared to the composite spectrum. We 
use LBQS, HST and SDSS QSO composite spectra given 
in Francis et al. (1991), Zheng et al. (1997) and Vanden 
Berk et al. (2001) respectively. First we shift the composite 
spectrum to the QSO emission redshift. For each extinc- 
tion curve we predict the reddened QSO spectrum by mul- 
tiplying the shifted composite spectrum by exp[— r(A)]. We 
mask the wavelength ranges of strong emission and absorp- 
tion lines. By varying the flux normalization and Av„, (or 
E(B-V)) we match the observed spectrum with our model 
reddened spectrum (see Fig. ^ . Results of our best fit mod- 
els with different extinction curves and associated reduced 
are summarised in Table U The error in Av^ includes the 
effect of errors in the parameters of the extinction curve. For 
both the systems the lowest value of is obtained for the 
dust extinction models with the extinction curve of LMC2- 
Supershell irrespective of our choice of composite spectrum. 



3 DISCUSSION ON INDIVIDUAL SYSTEMS 

3.1 The Zabs = 1.3265 absorption system towards 
J085042. 21+515911. 7 

The spectrum of this QSO (with u-K ^ 4.8 mag) shows a 
curvature around 2175 A in the rest frame of the Mg ll 
system (see top panel of Fig. [2}. The Milky Way extinc- 
tion curve that fits the 2175 A feature over predicts the 
QSO flux below 4500 A. This is also the case (not shown 
in Fig. [2)| when we use the average LMC extinction curve. 
The SMC extinction curve that flts the observed spectrum 
around 5500 A under predicts flux below 4500 A. The LMC2 
extinction curve that has a shallow UV-bump and relatively 
high extinction at rest wavelength Ar < 2100 A compared 
to that of the Milky Way (and of the LMC) reproduces the 
data better. We also fitted the spectra of ~ 200 non-BAL 
QSOs with Zem =1.892±0.005 using SMC extinction curve 
and find the probability of the 2175 A feature being pro- 
duced by QSO-to-QSO spectral variation to be <0.03. For 

1 details of MPFIT IDL routine can be found from 
|http://cow.physics.wisc.edu/~craigm/idl/titting.html| 
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Figure 2. The SDSS spectrum of J0850+5159 (top) and 
J0852+3435(bottom) are fitted with SDSS composite spectrum 
corrected using average extinction curves from the Milky Way 
(dotted), the LMC2 Supershell (solid) and the SMC (dashed). 
Rest wavelengths at z^bsa-re indicated at the top of the figure. 
The arrow marks the location of 2175 A feature at Zabs- For pre- 
sentation purpose the observed spectrum is boxcor smoothed by 
10 pixels. 



the range of composite spectra considered here, the mean 
Av, is 0.74±0.07 (and E(B-V) = 0.27). Our best fit model 
prediction J = 16.7, H — 15.8 and K — 14.7 mag, with 
a typical uncertainty of 0.2 mag, agrees well with the ob- 
served values, J = 16.99±0.23, H = 15.72±0.16 and K = 
15.27±0.15, from Strutskie et al. (2006). The excess flux 
predicted in K band can be mainly attributed to the known 
host galaxy contribution (^^0.3 mag) in the SDSS composite 
spectrum (See Vanden Berk et al. 2001) at Ar > 6OOOA. 
From Table 2 of Gordon et al. (2003) we have, 

iV(HI) ~ ^(6-97 ± 0.67) x lO^^cm"^ (3) 

Here k is the ratio of dust-to-gas ratio in the absorption sys- 
tem to that in LMC. We can use the equivalent widths of 
weak transitions (Table [Sjl to estimate the column densities 



of species and derive the depletion factors onto dust-grains 
assuming Zn is not depleted: —0.5,-0.9, —1.0, —0.9, and 
-1.0 for [Si/Zn], [Cr/Zn], [Fe/Zn], [Ti/Zn], and [Mn/Zn] 
respectively. Note that Zn llA2062 is blended with a Cr 11 
absorption but we can remove the contribution from Cr 11 
by using the average N{Cr 11) derived from unblended lines. 
Depletion is higher than what is typically seen in high- 2; ab- 
sorption systems and are intermediate between what is seen 
in the Cold and Warm phases in the Galactic ISM (Welty 
et al. 1999). Clearly more than 90% of the refractory metals 
are in dust as in the LMC. If we assume k ~ 0.9 then we 
get NiU i) = (5.73±1.10)xl0^^ cm'^ 

Our GMRT spectrum shows very strong absorption that 
can be modelled with two Gaussian components (see Table[2] 
and upper panel in Fig.[T}. The radio spectrum of this source 
is flat and VLBA maps at 5 and 15 GHz show that more 
than 90% of the flux density is in the unresolved core (Taylor 
et al. 2005). Using fc = 0.9 in Eq. 1, we derive iV(H l) = 
3.02±0.84 X 10" Ts cm"^ From the total N{R l) inferred 
from extinction we can estimate the spin temperature, Ts ~ 
190^gg* K. Note that this value is an upper limit as the 
dust-to-gas ratio may be larger than 1. All this is consistent 
with the system being associated with a cold neutral medium 
with dust properties similar to that seen in LMC2 supershell 
sample of Gordon et al. (2003). 

3.2 The Zs.hs = 1.3095 absorption system towards 
J085244.74+343540.5 

This QSO has u-K~5.6 mag and shows a curvature around 
2175 A in the rest frame of the Mg 11 system (see bot- 
tom panel of Fig. [2]) that is stronger than for the previ- 
ous quasar. The Mg 11 system at Zahs~ 1.3095 shows ab- 
sorption lines that are also stronger than that of the previ- 
ous system. The SDSS spectrum of this quasar is also well 
reproduced when applying a LMC2 like extinction curve 
to the QSO composite spectrum (with Ay^ = 1.00±0.09 
and E(B-V) = 0.36). Using spectra of ~ 330 non-BAL 
QSOs with Zem =1.655±0.005, we flnd the probability of 
the 2175 A feature being produced by QSO-to-QSO spec- 
tral variation to be <0.001. Our best fit model prediction J 

— 16.5, H — 15.2 and K — 14.2 mag agrees reasonably with 
the observed values, J = 16.75±0.14, H = 15.57±0.11 and K 

— 15.10±0.12, considering the possible uncertainties in the 
QSO composite spectrum discussed above. We derive lower 
limits on the column densities using the weaker metal lines. 
The absence of Cr llA2056 and Fe llA2249 features together 
with the high equivalent width of the Zn ll-|-Cr 11 blend at 
Xr = 2062 A is consistent with higher depletion factor and 
higher metallicity in this system. The detection of Ca 11 ab- 
sorption strongly supports this conclusion. From the SED 
fit and assuming k = 1 we derive N{R l) = 6.97±1.30 x 10^^ 
cm~^. This may well be an upper limit as k could be larger 
than 1. 

This system, despite having A'^(H i)k similar to that 
in the Zabs= 1.3265 towards J0850-I-5159, has a factor of 
two lower total integrated A'^(H l) from 21-cm absorption. 
This source has a fiat radio spectrum and is unresolved 
in the VLA A-array 8.4 GHz imaged having a resolution 

2 from CLASS (Cosmic Lens All-Sky Survey) and JVAS (Jo- 
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of 0". 255x0". 236. However, high resolution (few mas scale) 
VLBA map is not available for this source. Therefore, while 
/c could be close to 1 its exact value is still uncertain. We 
get a constraint iV(H l)/c/Ts = 1.31 ± 0.26 x 10^^ cm'^ 
K~^. By comparing the two A'^(H l) estimates we derive 
Ta/fc = 5361 gg* K. As fc can be less than unity and k is 
probably larger than one, the above value is an upper limit 
for the gas kinetic temperature. This is consistent with the 
gas being part of a cold neutral medium. 



4 DISCUSSION 

During our recently completed GMRT survey to search for 
21-cm absorption from intermediate redshift (1.10< z < 
1.45) Mg II systems, we discovered two red QSOs with strong 
z ^ 1.3 Mg II absorption systems along the line of sight. The 
observed spectral energy distribution of these QSOs are con- 
sistent with the QSO spectrum being reddened by dust in 
the intervening Mg ll systems. We fitted the SDSS spectrum 
using three QSO SDSS composite spectra reddened by dif- 
ferent extinction curves. We find that the dust properties 
of the absorbers are similar to what is seen in the LMC2 
supershell. In particular, we detect the 2175 A UV extinc- 
tion bump in both individual spectra. The neutral hydrogen 
column densities, A'^(H l), inferred from the extinction are 
consistent with the absorbing gas being of high column den- 
sity (i.e log N(H l)~21.7). We used weak metal absorption 
lines to estimate the column densities of ions of different 
species. Inferred metallicities are consistent with near solar 
values and depletion factors are similar to what is measured 
in the Cold Neutral Medium of LMC. 

These two quasars are rather weak (~ 50 mjy) in ra- 
dio emission and are not ideally suited for 21-cm absorption 
observations. We nonetheless detect 21-cm absorption from 
both systems thanks to the favorable physical conditions 
in the absorbing gas. These are the first detections of 21-cm 
absorption at high-z towards such faint background sources. 
The inferred spin temperatures in these systems are consis- 
tent with that of a cold neutral medium gas. 

The metal line equivalent widths and E(B-V) about 0.3 
measured in the two systems discussed here are higher than 
those derived from the SDSS composite spectra of Mg ll and 
Na I absorption systems (see York et al. 2006a; Wild, Hewett 
& Pettini 2006). The E(B-V) values are also higher than the 
median E(B-V) found for star forming galaxies at z = 2—3 
and a factor of 2 less than that measured in submillimeter- 
selected galaxies (see Figure 15 of Wild, Hewett & Pettini 
2006). The dust content and A''(H l) in the two systems dis- 
cussed in this paper are comparable to that observed in the 
21-cm system toward AO 0234-1-164. Thus these two systems 
are ideally suited for high resolution spectroscopic investi- 
gation of physical conditions in the interstellar medium of 
the corresponding absorbing galaxies. 
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